Control of spin in quantum dots with non- Fermi liquid correlations 
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Spin effects in the transport properties of a quantum dot with spin-charge separation are inves- 
tigated. It is found that the non-linear transport spectra are dominated by spin dynamics. Strong 
spin polarization effects are observed in a magnetic field. They can be controlled by varying gate 
and bias voltages. Complete polarization is stable against interactions. When polarization is not 
complete it is power-law enhanced by non-Fermi liquid effects. 

PACS numbers: 73.63.Kv, 71.10.Pm, 72.25.-b 



Spin phenomena in transport properties of low- 
dimensional quantum systems have become a subject of 
increasing interest QJ^ . Several fundamental effects have 
been predicted when controlling transport of electrons 
one by one in quantum dots as, for instance, spin block- 
ade due to selection rules [^ and parity effects in the 
Coulomb blockade |^,||]. There are also perspectives of 
applications in spin-electronics, quantum computing and 
communication j^. Previous works have been focusing 
on spin transport in two dimensional (2D) quantum dots 
connected to non-interacting leads and in the presence of 
a magnetic field |^] , including also an oscillating magnetic 
electron spin resonance component ||] . Spin transport in 
circuits with ferromagnetic elements and in the presence 
of a Luttinger-liquid interaction P-pj|| have been con- 
sidered. Of fundamental interest is the spin control of 
electron transport in the presence of correlations since in 
nanoscale devices the latter are very important. 

In the present paper, we derive a general theory for 
spin and charge transport through a quantum dot formed 
in a Luttinger liquid. We consider spin effects in the 
presence of a magnetic field. Specifically, we investigate 
to what extend non-Fermi liquid behaviour influences 
spin polarization. We find that spin-charge separation 
strongly affects the current- volt age characteristics. The 
spin leads to rich structure in the non-linear differential 
conductance that reflects both the collective spin density 
excitations and the orientations of the total spin in the 
quantum dot. A magnetic field in the quantum dot can 
spin-polarize the current strongly. This can be controlled 
by varying gate- or bias- volt ages. Full spin polarization 
can be achieved. Non-complete polarization is power-law 
enhanced by the non- Fermi liquid correlations. 

We start from a clean Luttinger liquid with spin. The 
charge interaction parameter is ^p = (l+2V(0)/7r'UF)~^^^ 
with vf the Fermi velocity and V{q) the Fourier trans- 
form of the electron interaction {h = 1). For the ex- 
change interaction, we assume g^ = 1. The inter- 
acting ID electrons are mapped, via bosonization, to 
a harmonic Hamiltonian |lj]. The low-energy excita- 



tions are charge and spin density waves, with disper- 
sions u;u{q) = VF\q\/gu = Vu\q\- Here, u = p^cr label 
charge (p) and spin (a). Spin-charge separation implies 
Vp 7^ Vfj. The slowly varying (on the scale of 27r/cp ) 
parts of the densities are given in terms of field operators 
6^(x), u{x) = p^{x) +p^p|(x) ^ uo ^ ^J2pii^J:Qy{x), 
with pp = -\-^ Pa = — and the mean charge and spin 
densities po = 2/cF/7r and ctq = 0, respectively. 

The quantum dot is formed by barriers {Ut/po)S{x—Xi) 
at positions xi < X2 given by the tunneling Hamiltonian 

^t Yl l[cos^{N-^iyod-a)cos^{N^-a) (1) 

q;=0,1 u 



with N^ = [6^(x2) ± e^(xi)]y277r. Physically, N' is 
the deviation of the number of electrons from the mean 
value, no = dpo^ in the interval d = X2 — xi. The excess 
charge then is Q = —eN~, and the component of the 
total spin of the electrons parallel to the quantization 
axis (assumed parallel to the ID system) S = N~ /2. The 
numbers of imbalanced electrons and spins between the 
leads are N^ and A/'+, respectively. The coupling to the 
source-drain bias V and the gate voltage Vg is described 
by Hv = -e(V7V+/2 + l/g7V-(5), with S the ratio between 
gate and total capacitances. The effect of an external 
magnetic field is described by a local Zeeman term in 
the region between the barriers, Hb = —gBPBBN~ /2 = 
-EbN~/2, with the Lande- factor gB 0,§. 

The currents are calculated as the stationary lim- 
its of transferred charges and spins I^ = I^ -\- p^Ii = 
(e/2) limt^oo(A^j^(^))- The brackets include both ther- 
mal and statistical averages over the collective modes at 
X 7^ xi,X2 with the density matrix reduced to N^ [|[. 

For obtaining the non-linear current- voltage character- 
istics we consider the dynamics of the system described 
by the variables N^ under the influence of the external 
fields in the 4D periodic potential Eq. (^. For high 
barriers, tunneling between nearest-neighbored minima 
dominate, with amplitude A that is related to Ut via 
the WKB-approximation pj. These correspond to pro- 



cesses A^^ 



N- ± 1 and N- 



N^ zb 1 associated 



with changes of charge and spin numbers in the dot, re- 
spectively, and A^^ -^ N^ ± 1 that transfer current. We 
consider sequential tunneling with the temperature much 
smaller than the dot level spacing. This can be described 
by a master equation for charges and spins with rates. 
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EiE)= y^WP{ep)W:^{e,Ws\E-nep-me,) , (2) 



where the tunneling rate of a single barrier 



v(9) 



<-)^£;(^) 



1-1/ 



r(i/5) 



i(3E\ 2 



(3) 



depends on the effective interaction parameter g/2 = 
9p9a/{9p-^ 9a) and the frequency cutoff a;c. The weights, 
W^(ejy), at the discrete energies e^^ are {P~^ = k^T <^ e^) 



^n(^.) 



1/2^. 



n\T{l/2g,) 



-e.n/a;e5,(^) . (4) 



In order to understand the rather complex behaviour 
of the transport spectra it is useful to recall the charac- 
teristic energy scales in (|2|). These are the discretization 
energies corresponding to charge and spin modes in the 
quantum dot relative to the energy of the ground state. 



e^ =ujj,{q = n/d) = 2g^E^ , 



(5) 



with the addition energies for charge and spin 

E, = 7rv^/2gld. (6) 

Without interaction, the addition energies are Ep = 
Ea^ = 7rvY/2d = £^p 7^ 0, due to the Pauli principle, and 
the discreteness of the dot levels. On the other hand, for 
strong Coulomb interaction, Ep oc V{q ^ 0) ;» £^cr ||]. 

The addition energies determine the ground state en- 
ergy of n charges with the half-integer total spin S = 
Sn/2, Eo{n) = [Ep (n - n^f + E^sl - EBSn]/2. The 
reference particle number n^ = eVgd/Ep + no is defined 
by the gate voltage. The energy differences of the many- 
body states of n + 1 and n electrons are 

E W 

fi{n, s, /, m) = ^ (2n + 1 - 2ng) + ^(4+i - 4) 

Eb 

^(^n+l - Sn) + lep + mCa • (7) 

Positive or negative integers / and m denote the differ- 
ences of the numbers of charge and spin excitation quanta 
with the energies (||), respectively. These do not change 
neither the number of particles nor the total spin in the 
quantum dot. The energy differences fi{n,s^l,m) play 
the role of chemical potentials of the dot and define the 
transport regions. For symmetric bias, for instance, the 
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FIG. 1. Differential conductance dl/dV as a function of 
the source-drain bias eV/Ep for gp = 0.4 (cjc = lO^^cr, units 
10-^(cJc/^cx)^^"^^/^^e^AV4cj?). Top: Eb = 0, rig = 0.58; 
bottom: Eb = 0.4 ^cr, rig = 0.548. 



condition V/2 > jiin^ s^l^m) > —V/2 defines the allowed 
transport channels. For V ^ and Eb <C Ecr one finds 
the Coulomb blockade peaks at gate voltages Vg that 
correspond to ii{n^s^l^m). However, due to the spin, 
the separation of the peaks depend on the parity of n, 
/^yn+i,n ^ [Ep^{-lY+^{E^-EB)]/e5. Without inter- 
action, one has to replace here Ep = E^ = Ep, in order 
to get the separation of the linear conductance peaks. 

In the following, we consider the limit T = 0. Re- 
sults for the differential conductances d//dV for Eb = 
and Eb = OAE^ as functions of V are shown in Fig. [^. 
Zero bias voltage has been assumed at the position of a 
conductance peak corresponding to an n(even)-to-(n + l) 
ground-state-to-ground-state transition. The differential 
conductance shows sharp peaks at bias voltages Vnism 
at which a new transport channel enters the above bias 
voltage window. Above Vnism^ the conductance as a func- 
tion of V drops according to the interaction-induced non- 
Fermi liquid power law (V — VnismY^^~'^ - 

From (0) one can easily identify the spectral origins 
of the peaks in the conductance spectra. At low volt- 
ages, the spectra are completely dominated by spin ex- 
citations due to spin-charge separation. The discretiza- 
tion and addition energies corresponding to the spin are 
factors gp and ^^, respectively, smaller than those corre- 
sponding to the charge. For 5 = peak (e) corresponds 
to a charge density excitation at e^, while (/) is due to 
the ground-state to ground-state transition at Ep — E^. 
All of the other peaks in Fig. |l| are spin-related. Be- 
cause ga = I the transition at 2Ea is degenerate with 
the spin density excitation at Ca- (peak (a), with multi- 
ples (d) and {g)). A finite exchange would remove this 
degeneracy and discriminate between spin addition ener- 
gies and spin density waves. The two small features {b) 
and (c) are combinations of the excitations (e) and (/) 
with E(j, they corresponds to Cp — 2E(j and Ep — 3£^cr, re- 
spectively. In the non-interacting limit, the peaks in the 
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FIG. 2. Density plot of the differential conductance in the 
(ey/^p,ng)-plane for interaction gp = 0.4 and magnetic field 
Eb = OAEa (uJc = lO^^cT, color code (left) with units 
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non-linear differential conductance appear at bias volt- 
ages that are multiples of AEp due to the absence of spin 
charge separation. Apart from the peak at F = 0, all of 
the peaks are absent without electron interaction, for the 
bias voltages in Fig. ^ (top). 

The spin-related features are even more strikingly dis- 
played in the spectra for Eb = OAE^ (Fig. |l| bottom). 
All of the peaks in Fig. ffl (top) acquire Zeeman side bands 
corresponding to energies ^peaki^B- Exceptions are {b) 
and (c). They only have sidebands £^peak + Eb since the 
initial states corresponding to the lower sidebands cannot 
be occupied by electrons entering the dot at T = 0. 

As an example of the rich spin-related structure in the 
transport spectrum a density plot of the differential con- 
ductance in the {eV/Ep^ng)-p\dine for Eb = OAEa- is 
shown in Fig. g. The black regions correspond to the 
Coulomb blockade. When increasing the bias voltage the 
the excitation spectrum displays a considerable number 
of spin-related transitions. With non- vanishing exchange 
interaction, ^^r 7^ 1, accidental degeneracies due to spin 
addition and excitation are lifted and the spectrum be- 
comes even more complex. 

Figure H shows the behaviour of a current peak for 
fixed bias, /(%), when changing the magnetic field. As 
B is changed, peak height and position vary with periods 
2Ecr and 4£^cr, respectively (Fig. || bottom). This can be 
understood by considering the processes that contribute 
towards the current. We start by discussing the peak 
position for F ^ 0. For small B and keeping rig as to 
match the maximum of the peak (Fig. top) one finds 
from (^ that with increasing B one has to adjust Ug 
to lower values ex —Eb/'^ since 5^+1 — 5^ = +1 which 
corresponds to the s^ = ^ ^5^+1 = 1 transition. When 
Eb > 2^cr, the energy of the state s^ = 2 becomes lower 
than that of the 5^ = state such that transport gets 
now support from transitions 5^ = 2 ^ 5^+1 = 1 with 
Sn+i — Sn = —^ while 5n+i + 5^ = 3 instead of 1. When 
increasing B further, rig has now to be adjusted to higher 
values (X -^EbJ'^ in order to compensate for the Zeeman 
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FIG. 3. Top: Current /, normalized to the value at B = 0, 
as a function of the gate voltage rig for different magnetic 
fields Eb/E^ (inset) for gp = 0.4 and at eV = 0.1 ^p 
(cjc = 10^ ^cj)- Bottom: position of maximum of peak, no 
(dashed, right scale) and current at peak maximum, /o (full 
line, left scale) as a function of Eb/ Ea. 



shift. The original position is reached after a total change 
of B corresponding to Eb = ^E^. For ]/ ^ 0, the peak 
position remains unchanged as long as energies of the 
transitions are inside the interval (— V/2, V/2). 

For understanding the peak height we consider V 7^ 0. 
The bias V gives the width of the current peak. First, 
we remember that for Eb <C eF, the current is due to 
two transitions, namely 5^ = ^ s^+i = ±1. This 
leads to the asymmetry of the (non-spin polarized) peak 
for Eb = that can be observed in Fig. || (top). When 
Eb ^ eV the contribution of the ground-to-excited-state 
transition s^ = ^ 5^+1 = — 1 is suppressed. Further 
increasing B, the current peak becomes symmetric, com- 
pletely polarized (cf. Fig. Q), its height is reduced and 
remains constant. For Ecr <C Eb < 2Ecr transport gets 
support also from transition 5^ = 2 ^ s^+i = 1 such 
that the peak height starts to become again asymmetric 
and to increase until, exactly at Eb = ^E^j, both contri- 
butions are equally important. Then, the current peak 
acquires the same shape as for Eb = 0, but reflected at 
rig = 1/2. Increasing B produces oscillatory behavior 
due to further changes of the values of Sn and 5^+1- The 
current has its maximum value at rig = 1/2 whenever 
Eb is an odd multiple of E^^. This fully polarized cur- 
rent states are also reflected in the digital behaviour jl^] 
displayed in Fig. ^ (bottom). The latter becomes less 
stable when increasing the bias voltage V. 

Figure El shows the behaviour of the current polariza- 
tion, P = la/Ip' The top- left panel shows P as a function 
the source-drain voltage for a given magnetic field and 
different interaction strengths, with rig at the maximum 
of the linear conductance. The polarization is complete 
for eV < 2Eb^ independently of the interaction. When 
eV > 2Eb, P decreases as a function of V according to 
a non-Fermi liquid power law and is higher for stronger 
interaction. Thus, correlations enhance the polarization 
when the latter is not complete. 
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FIG. 4. Current polarization P 



for different inter- 



actions Qp (inset) as a function of cV/Eb for Eb — ^AE^ 
and rig at the maximum of the conductance peak (top left); 
P in the plane {eV/Ep^Ug) for Qp — 0.4 and magnetic field 
Eb — 0.5 ^cr (top right), Eb — Ea (bottom left) and 
Eb = 1.5 E(j (bottom right, includes color code). 



tures, one would expect that coherent tunneling processes 
dominate. Including spin, these lead to the well-known 
quantum-dot Kondo physics which is not considered in 
our approach [Q. Therefore, the above results apply to 
temperatures higher than the Kondo temperature. 

Our results for the transport spectra are consistent 
with several of the non-linear features observed recently 
in a one dimensional (ID) quantum dot formed by two 
impurities in a cleaved-edge overgrowth quantum wire 
[|l6[. More detailed experiments, however, are needed in 
order to test the above predictions, especially concern- 
ing the control of the spin. We expect that the effects 
can be observed in the transport through double barriers 
formed in cleaved-edge-overgrowth quantum wires, and 
in carbon nanotubes p/7|] . 

Acknowledgements: This work has been supported by 
Itahan MURST via PRIN 2000, by the EU within TMR 
and RTN programmes, and by the DFG. 



The other three panels of Fig. Q show density-plots of 
the polarization for fixed interaction. Details of the be- 
haviour can be understood from the discussion related to 
Fig. ||. Varying the magnetic field, complete polarization, 
P = 1 (yellow), is transferred between different regions 
of the (eV, ng)-plane. The top-right panel corresponds to 
Eb = 0.5 E (J were the dominating Sn = ^ s^+i = 1 
channel leads to complete spin polarization near the 
left-hand edge of the region of non-zero current. Near 
Eb = Ecr (bottom-left) the current peak is symmetrically 
spin polarized (for small F, yellow diamond). Increasing 
the magnetic field fully polarizes the right-hand edge of 
the current region {Eb = 1.5£^cr, bottom-right). Exactly 
at Eb = 2Ea- the two involved transport channels have 
equal weight and spin polarization is exactly zero, as for 
Eb = 0. By further increasing B, the behaviour of the 
polarization is reversed, P = — 1. It can be displayed by 
the same panels, but in opposite direction. The period- 
icity of P with respect to Eb corresponds to AE^^. 

In conclusion, we have shown how one can control the 
spin properties of the transport through a ID quantum 
dot embedded in a non-Fermi liquid by changing mag- 
netic field, and gate- and source-drain voltages. This is 
due to the interaction that separates the energy scales of 
the charge and the spin excitations such that E^ <^ Ep. 
Complete spin polarization can be achieved in spite of 
the presence of correlations. Once it is achieved, it is not 
influenced by the interaction. When polarization is not 
complete, it is enhanced by the non- Fermi liquid corre- 
lations. This shows that the electron spin is crucial for 
understanding non-linear transport in ID quantum dots. 

The above results have been obtained for T = 0. When 
T > 0, we expect temperature-induced Luttinger liquid 
power-law broadening of the conductance peaks, and cor- 
respondingly a smearing of the spin polarization features 
which is the subject of future work. At very low tempera- 



[1] G.A. Prinz, Physics Today 282, 58 (1995); Science 282, 

1660 (1998). 
[2] J.M. Kikkawa, and D.D. Awschalom, Phys. Rev. Lett. 
80, 4313 (1998); D.D. Awschalom, and J.M. Kikkawa, 
Phys. Today 52, 33 (1999). 
[3] D. Weinmann, W. Hausler, and B. Kramer, Phys. Rev. 

Lett. 74, 984 (1995); 
[4] S. Tarucha, D.G. Austing, T. Honda, R.J. van der Hage, 
and L.P Kouwenhoven, Phys. Rev. Lett. 77, 3613 (1996). 
[5] T. Kleimann, M. Sassetti, B. Kramer, and A. Yacoby, 

Phys. Rev. B 62, 8144 (2000). 
[6] A. Imamoglu, D.D. Awschalom, G. Burkard, D.P. DiVin- 
cenzo, D. Loss, M. Sherwin, and A. Small, Phys. Rev. 
Lett. 83, 4204 (1999). 
[7] P. Recher, E.V. Sukhorukov, and D. Loss, Phys. Rev. 

Lett. 85, 1962 (2000). 
[8] H.A. Engel, and D. Loss, Phys. Rev. Lett. 86, 4648 

(2001). 
[9] A. Brathaas, Y. Nazarov, and G.E.W. Bauer, Phys. Rev. 
Lett. 84, 2481 (2000). 
[10] Q. Si, Phys. Rev. Lett. 81, 3191 (1998). 
[11] L. Balents, and R. Egger, Phys. Rev. Lett. 85, 3464 

(2000). 
[12] J. Voit, Rep. Progr. Phys. 58, 977 (1995). 
[13] A. Braggio, M. Grifoni, M. Sassetti, and F. NapoU, Eu- 

rophys. Lett. 50, 236 (2000). 
[14] A. Furusaki, Phys. Rev. B 57, 7141 (1998). 
[15] M. Ciorga, A.S. Sachrajda, P. Hawrylak, C. Gould, P. 
Zawadzki, S. Jullian, Y. Feng, and Z. Wasilewski, Phys. 
Rev. B 61, R16315 (2000). 
[16] O.M. Auslaender, A. Yacoby, R. de Picciotto, K.W. Bald- 
win, L.N. Pfeiffer, and K.W. West, Phys. Rev. Lett. 84, 
1764 (2000). 
[17] H.W.Ch. Postma, M. de Jonge, Z. Yao, and C. Dekker, 
Phys. Rev. B 62, R10653 (2000). 



